Testing the meson cloud in the nucleon in Drell-Yan processes 
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Abstract 

We discuss the present status of the u — d asymmetry in the nucleon and 
analize the quantities which are best suited to verify the asymmetry. We find 



^ ■ that the Drell-Yan asymmetry is the quantity insensitive to the valence quark 



distributions and very sensitive to the flavour asymmetry of the sea. We com- 
pare the prediction of the meson cloud model with different experimental data 
including the Fermilab E772 data and recent data of the NA51 Collaboration 

at CERN and make predictions for the planned Drell-Yan experiments. 
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I. INTRODUCTION 



The deviation of the Gottfried Sum Rule from its classical value [|T| 

Sg = fmx)-F^x)]-= l - (1) 

JO X 6 

observed by the New Muon Collaboration (NMC) at CERN 0,|3| has created a large interest 
on the possible sources of the violation. In the NMC experiment the neutron structure 
function which enters the sum rule is deduced from deep inelastic scattering off deuterium. 
It could be biased by nuclear two-body effects which were ignored in the NMC analysis. 
While shadowing effects 0-0 cause the real value of the Gottfried Sum Rule to be even 
smaller than the value given by NMC, the anti-shadowing effect, due to the presence of 
virtual mesons which bound the deuteron, tends to restore the classical value |||7[]. A 
recent unfolding of the shadowing and anti-shadowing effects |7[ suggest, however, a genuine 
violation. This can be understood as consequence of an internal asymmetry d > u of the 
quarks in the proton (the opposite asymmetry is expected for the neutron if the proton- 
neutron charge symmetry holds). The asymmetry has been confirmed recently by the NA51 
collaboration group at CERN ||. Since at large Q 2 the perturbative QCD evolution is 
flavour independent and, to leading order in logQ 2 , generates equal number of uu and dd 
sea quarks nonperturbative effects must play an important role here. The effect of the 
asymmetry has been predicted by the models in which the physical nucleon contains an 
admixture of the irN and 7rA components in the Fock expansion |10j. The predicted effect 



of the asymmetry agrees [TT|] with that deduced from the NA51 experiment. 

The idea of the u — d asymmetry is not new. It was considered a decade ago by Ito 
et al. ]TJ] as a possible explanation for the slope of the rapidity distribution of dilepton 
production in the proton-Pt collision at Fermilab. An alternative interpretation stimulated 
by the discovery of the EMC effect invoked the enhancement of the nuclear sea in Pt with 



respect to a collection of free nucleons ||13|| . In general both effects can coexist 



In the advent of high precision data on the dilepton production in the proton-proton and 



proton-deuteron scattering [[[J], it is interesting to review the present status of our knowlegde 



on the u — d asymmetry. We will discuss the quantities which should most unambiguously 
confirm the asymmetry and allow for verification of different theoretical concepts. We also 
confront the prediction of the meson cloud model with the existing data. 
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II. GOTTFRIED SUM RULE 



The Gottfried Sum Rule M (GSR) addresses the value of the integral over x of the 
difference of the F 2 (x) structure function of the proton (p) and neutron in). It is written^] 

as 



1 [Fi(x) - F?(x)] - 



x 



1 4 



o 9 
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<0) - <0) 



d;(x) - di{x) 



i 

+ 9 

+ 2 \ ^ [u p {x) - u n (x)} - - dp(x) - d n (x) 



(2) 



dx 



where u v p (x) 
further to 



u p (x) — u p (x), etc. Baryon number conservation reduces the expression 



s = i + " 



7 (x) - u n (x)\ - - d p (x) - d n (x) 



dx . 



(3) 



As seen from Eq.(|3|) the valence quarks do not influence the Gottfried Sum Rule violation 
although they are of crucial importance for the GSR integrand. Assuming further charge 
symmetry for the nucleon ) = d n (x), etc. and making the customary assump- 

tion that u p (x) = d p (x), one finds the classical value of 1/3. The NMC experiment [0] of the 
relevant structure functions, over the interval 0.004 < x < 0.8, yielded when extrapolated 
to < x < 1, 



f 1 d,T, 
/ [F%(x) - F^{n)} — = 0.24 ±0.016 , 
J oc 



(4) 



at Q 2 = 5 GeV 2 . It should be noted that QCD corrections do not play any role here. While 
the leading order corrections to the Gottfried Sum Rule cancel, the higher order corrections 
are negligibly small. 

The integrand of the Gottfried Sum Rule has been obtained from the ratio R = F£ / F% 



and the deuteron F^ structure function 



x The structure functions F2(x) and the quark distribution functions q(x) are, of course, func- 
tions of Q 2 . However, the Q 2 dependence is suppressed to keep the expressions from being too 
cumbersome. 



Ff(x)-F 2 >)=2F 2 d (x)[ T |||. (5) 

In the first evaluation of GSR the deuteron structure function was taken from a global fit to 
the results of earlier experiments. In the meantime NMC has published values of F% and F 2 d 



of its own accurate measurements at low x (T^|. The old value of GSR M has been updated 
ito 

\ -F 2 x (n)] — = 0.235 ±0.026 . (6) 

The total error here is larger than in (^) due to a more extensive examination of the sys- 
tematic uncertainties. The quoted error bar does not include the effects of shadowing and 
antishadowing. 

In the most general case not only the so-called SU(2)q |16| charge-symmetry is vio- 
lated but also the isospin symmetry between proton and neutron SU(2)i. In this case the 
antiquark distributions in proton and neutron can be expressed as: 

u p (x) = q{x) - ^Aq(x) - -A/(x) , 
d p {x) = q(x) + ^A Q (x) - -A/0) , 
u n (x) = q(x) + ^Aq(x) + ^A 7 (x) , 

d n {x) = q{x) - ^A Q (x) + -Ajix) . (7) 

The signs in front of Aq(x) and A/(x) have been chosen to assure positivity of Aq and 
Aj in the case the asymmetries give the effect required by the NMC result. With the 
parameterization (^) the Gottfried Sum Rule @ can be written as 

The factor 10/9 shows the sensitivity to the SU(2)j symmetry violation. The violation 
of the Gottfried Sum Rule exclusively via breaking the SU(2)q symmetry requires Aq = 
0.14 ± 0.03. On the other hand violation of GSR exclusively via breaking the SU(2)j 
symmetry requires Aj = 0.08 ± 0.02. 



It has been argued in Ref. |T7| that both the effects of asymmetric £77(2) sea and the 
proton-neutron isospin symmetry breaking will be very difficult to disentangle as they may, 
in general, lead to very much the same behaviour both in deep inelastic scattering and Drell- 
Yan processes. A careful analysis [IIJ based on the cr-term and scale Ward identity suggests 
rather Aq ^> A/. Therefore in the further analysis we will neglect the effect of the ££7/(2) 
symmetry violation. 



Since the standard Altarelli-Parisi |18[ evolution equations generate equal number of 
uu and dd pairs, one does not expect strong scale dependence of the GSR. The two-loop 
evolution gives a rather negligible effect The Pauli exclusion principle leads to some 
interference phenomena which produces only a small asymmetry ||19|| . 

The answer likely lies with more complicated nonperturbative physics. The presence of 
the pion cloud in the nucleon gives a natural explanation of the excess of d over u. It has 
been extensively analyzed in a series of papers fl2"0"|-f26|j with the result for the GSR being 
dependent on the details of the model, especially on vertex form factors. Restricting the 
form factors by fitting to the cross sections for high-energy production of neutrons and A ++ 
yields the GSR f2JJTg|27|] in rough agreement with that obtained by NMC 

The flavour content of the sea in the nucleons can be tested in Drell-Yan experiments, 
which is the subject of the present paper. In the first phenomenological analyses the anti- 
quark distributions were typically defined as 

u(x) = q(x) — ^A(x) and d(x) = q(x) + ^A(x) . (9) 

Different assumptions on A(x) lead to different predictions for the Drell-Yan production 
rate. For instance, the initial parameterization [^HJ used A(x) = A(l — x) k , which placed 



the d — u difference at large x (x > 0.05), led to very large values for d p (x)/u p (x) for x > 0.1 
These large ratios have been ruled out by a recent reanalysis of earlier Drell-Yan data p9 



More recent |30| parameterizations have a form A(a;) = A/s,x v (l — x)' ns (l + r fx), which places 
the bulk of the difference at smaller x. 

It is expected that future QCD lattice gauge calculations will be able to generate the 
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quark structure of the nucleon and evaluate the d p (x) and u p (x) distributions. At present one 
has to rely on phenomenological models consistent with our knowledge in other branches of 



hadronic physics. The meson cloud model [|3T|J27|1 seems to satisfy this criterion. Furthermore 
it is worth mentioning in this context that importance of pion loop effects for nucleon 
properties has been demonstrated in a recent lattice QCD calculation [^2] . 
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III. MESON CLOUD MODEL OF THE NUCLEON 



In this section we briefly sketch the meson cloud model (MCM) of the nucleon [26,31,27 
and present its prediction for the asymmetry of the light sea antiquarks. In this model the 
nucleon is viewed as a quark core, termed a bare nucleon, surrounded by the mesonic cloud. 
The nucleon wave function can be schematically written as a superposition of a few principle 
Fock components (only ttN and tcA are shown explicitly) 



\P)phys — V Z \p} CO re 

+ J dyd 2 k±(/) Nn (y,kA.)(]J^\pTT ;y,kA_) + \J^\nn + ;y, kj_)) 

+ jdy d 2 k ± <p Aw (y, k ± ) (y^l A++7T-; y, k ± ) - ^| A+vr ; y, k ± ) + ^| A 
+ ... 



~;y,k±) 



(10) 



with Z being the wave function renormalization constant which can be calculated by impos- 
ing the normalization condition (p\p) = 1. The <p(y, k±)'s are the light cone wave functions 
of the ttN, 7rA, etc. Fock states, y is the longitudinal momentum fraction of the 7r (meson) 
and k± its transverse momentum. 

It can be expected that the structure of the bare nucleon (core) is rather simple. Presum- 
ably, it can be described as a three quark system in the static limit. Of course, in the deep 
inelastic regime at higher Q 2 additional sea of perturbative nature is created unavoidably 
by the standard QCD evolution. 

The model includes all the mesons and baryons required in the description of the low 
energy nucleon-nucleon and hyperon-nucleon scattering, i.e. the 7T, K, p, ui, K* and the N, 
A, S, A and S*. In contrast to other approaches in the literature the model ensures charge 
conservation, baryon number and momentum sum rules by construction. 

The main ingredients of the model are the vertex coupling constants, the parton distribu- 
tion functions for the virtual mesons and baryons and the vertex form factors which account 
for the extended nature of the hadrons. The coupling constants are assumed to be related 



S 



via SU(3) symmetry which seems to be well established from low-energy hyperon-nucleon 
scattering. 

It was suggested in Ref. |?5] to use the light cone meson-baryon vertex form factor 

M 2 MB (y,ki)-m%~ 



F{y, k ± ) = exp 



?A 2 



(11) 



where k± is the transverse momentum of the meson and Mmb{u, k\) is the invariant mass 
of the intermediate two-body meson-baryon Fock state, 

MUv, kl) = 4^ + ■ (12) 

By construction, such form factors assure the momentum sum rule [|25|j3T[| . The parameters 
Amb are the principal nonperturbative parameters of the model. They have been determined 
from an analysis of the p — > n, A, A fragmentation spectra |27[ using light cone flux functions 
|27|] (A^ = 1.08 GeV 2 and A^ A = 0.98 GeV 2 ). With these parameters the pion exchange 
model gives a very good description of the spectra. 

In practice the probability of the Fock components with strange particles is rather small. 
For instance the probability to find a KA state in the nucleon is about 1%, whereas the 
probability to find a ttN state is 0.18. In all applications in the present paper the higher 
Fock states involving strange particles can be neglected. 

The x dependence of the structure functions in the meson cloud model can be written 
as a sum of components corresponding to the expansion given by Eq . (|10|) . 



F 2 N (x) = Z 



F 2 N core (x) + j:(5^F 2 (x)+5^F 2 (x 

MB 



(13) 



The contributions from the virtual mesons can be written as a convolution of the meson 
(baryon) structure functions and its longitudinal momentum distribution in the nucleon BB 



5^F 2 (x)= f 1 dyf M (y)F 2 M (-). (14) 
Jx y 

The same is true for the virtual baryons 

5 {B) F 2 (x)= f 1 dyf B (y)F 2 B core (*). (15) 
Jx y 
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In a natural way /m(z/) and are related via |3T 



f B (y) = f M (l - y) . (16) 

Eq. (0) (and also Eq. |T|) can be written in an equivalent form in terms of the quark 
distribution functions 

**QfM=[ l Mv)tf&-- ( 17 ) 
Jx y y 

The longitudinal momentum distributions (splitting functions, flux factors) of virtual mesons 
(or baryons) can be calculated assuming a model of the vertex and depend on the coupling 
constants and vertex form factors. Further details can be found in Refs. pipTfl . 

The parton distributions "measured" in pion-nucleus Drell-Yan processes are used 
for the mesons. The deep-inelastic structure functions of the bare baryons, F^ core (x,Q 2 ), 
Ff core (x,Q 2 ) are in principle unknown. In practical calculations it is usually assumed 
F 2,core( x iQ 2 ) = F 2,ph ys { x iQ 2 ) which is not full y consistent. Recently |36 



we 



have extracted F^ core by fitting the quark distributions in the bare nucleon, together with 
(parameter-free) mesonic corrections, to the experimental data on deep-inelastic scattering: 

(a) F|(x,Q 2 )-F 2 "(x,Q 2 ) §, 

(b) F?(x,Q 2 )/FZ(x,Q 2 ) i, 

(c) F 3 ^(x,g 2 ) n, 

(d) q(x,Q 2 ) @. 



The following simple parameterization has been used for the quark distributions in the bare 
proton at the initial scale Qq = 4 (GeV/c) 2 . Note, that we have used u-d symmetric sea 
quark distribution for the core. 

xu v , core (x, Ql) = N u x°- 38 (1 - x) 2A9 (l + 10.5x) , (18) 
xd v , core (x, Ql) = N d x om (l - x) 4 - 63 (l + 150x) , 
xS core (x,Q 2 ) =0.17(1 -x) 13 - 8 , 
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where 



^core ^s,core ^s,core d score ds,core 2,S SjCore 2s SjCore . (19) 



The details of the fit and a comprehensive discussion of DIS will be given in Ref. |36| . 
An example of the fit to F%(x) — F£{x) is shown in Fig.l. The resulting total sea quark 
distribution is compared in Fig. 2 with an experimental sea quark distribution obtained 



from (anti)neutrino induced reactions [p8|1 . The so-extracted parameterization of the quark 
distributions in the bare baryons can be used to calculate the cross-sections for both the 
lepton deep-inelastic scattering and for the Drell-Yan processes. In the present article we will 
present predictions(l) for the Drell-Yan dilepton production in elementary nucleon-nucleon 
collisions as well as for the nucleon-nucleus collisions. 
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IV. DRELL-YAN PROCESSES 



The Drell-Yan process [3£| involves the electromagnetic annihilation of a quark (an- 



tiquark) from the incident hadron A with an antiquark (quark) in the target hadron B. The 
resultant virtual photon materializes as a dilepton pair (£ + £~) with muons being the pair 
most readily detected in experiments. 

Ellis and Stirling have shown that the measurement of Drell-Yan cross sections in the 
proton-proton and proton-deuteron collisions provides information on the d p (x) /u p (x) ratio 



The cross section for the DY process can be written as 
da AB 4na 2 



EeJ [A^WbM+^M^M] , (20) 



dx\dx2 9sx\X2 j 
where s is the square of the center-of-mass energy and x± and x 2 are the longitudinal momen- 
tum fractions carried by the quarks of flavour /. The c[ A {x\) {<1a( x i)) an d Qb( x 2) (Qb{ x 2)) 
are the (anti-)quark distribution functions of the beam and target, respectively. The factor 
K(xi,x 2 ) accounts for the higher-order QCD corrections that enter the process. Its value 
over the kinematic range where experiments are carried out is typically 1.5. The values of 
X\ and x 2 are extracted from experiment via 

M 2 = sx x x 2 ~ 2P e +P e -{l - cos9 i+i -), (21) 

where M is the mass of the dilepton pair, P^+ and Pg- are the laboratory momenta of 
the leptons, and Qi+i- is the angle between their momenta vectors. The total longitudinal 
momentum of the lepton pair (P^+ + Pe-)i fixes Xi — x 2 via 

Xl -x 2 =x F = 2{Pt¥ + Pt - )L -l. (22) 

s 

In order to avoid spurious contributions to the DY yield from vector meson decays, all 
measurements are made for M > 4 GeV, and the region 9 < M < 11 GeV is excluded to 
avoid the T resonances. 
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The absolute value of the Drell-Yan cross section is biased by the uncertainty in ex- 
trapolating from time-like to space-like values of Q 2 when relating the Drell-Yan with deep- 
inelastic scattering which involves the factor K (see Eq.^y). In order to avoid the uncer- 



tainty it is desirable to consider ratios fli3| , |2l||T0|| rather than the absolute cross sections. 



Whether the i^-factors for the pp and pn collisions are identical can be checked by calcu- 
lating higher-order QCD corrections. In Fig. 3 we show the .fT-factors calculated according 
to the formalism developed in Ref. |34| for two different leading order quark distributions 
40|j4l| . Although the .fT-factor depends on the input quark distribution as well as on x\ 



and X2 its value is practically identical for proton-proton and proton-neutron collisions. The 
approximate equality K pp = K pn allows us to neglect the higher order corrections and limit 
ourselves to the much simpler leading order analysis. 

The quark distributions found from the procedure described in the former section can 



be further tested by comparison with the Drell-Yan E772 data |29| for the differential cross 
section M 3 d 2 a / 'dxpdM for the dilepton production in the p + d collision. By fitting the 
i^-factors the result of the calculation can be compared with the experimental data. In 
order to check the sensitivity to the u — d asymmetry we have also performed the calculation 
with symmetrized sea: 

u s (x) = u s (x) = d s (x) = d s (x) = . (23) 

In Fig. 4 we show the result of the fit for our original model (solid line, y 2 /N = 1.15) and the 
results obtained with the symmetrized (Eq. ^) sea distribution (dashed line, y 2 /N = 1.42). 
Altough there is some sensitivity to the d — u asymmetry, it can be easily compensated by 
a slightly different normalization factor. For comparision we also show the result obtained 
with the MSR(S' ) quark distributions (dash-dotted line, y^/N — 1-84). 

The present experimental data for the Drell-Yan processes in the elementary nucleon- 
nucleon collisions suffer from rather low statistics. Therefore at present one is forced to 
compare a theoretical calculation with the proton-nucleus experimental data. In the first 
approximation the cross section for the production of the dilepton pairs in the proton-nucleus 
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scattering can be expressed in terms of the elementary pp and pn processes as 

</ = Zo% + No% . (24) 



It has been shown [f|^,|TI]] that the ratio of the cross section for the scattering of protons from 
the nucleus with N — Z ^ to that from an isoscalar target such as deuterium is sensitive 
to the d p (x) —u p (x) difference. These ratios have been measured by the E772 Collaboration 
at FNAL |42| for carbon, calcium, iron and tungsten targets. Neglecting nuclear effects, 
elementary algebra leads to the following result for the ratio 

R = 2a DY (p + A) = 2Z N-Z 2a* n (x 1 ,x 2 ) 
DY Aa DY (p + d) A A aPP(x h x 2 ) +aP n (xi,x 2 ) 7 

where Z, N, A are number of protons, neutrons and the atomic number, respectively. In 



the large x 2 (target) limit the ratio takes a very simple form [29 



/ \ N-Z A(x) , . 

= 1+ a WfTWY (26) 

showing that the Drell-Yan processes with non-isoscalar targets are relevant for the issue of 
the asymmetry. 

The experimental ratios are consistent with symmetric quark distributions |2"P|, 



sec 



Fig. 5), which renders those data useless for establishing the asymmetry. Moreover, using 
asymmetric quarks distribution functions (solid and dashed lines) has a rather small effect 
on the ratio. The ratio obtained with the recent MSR(A) quark distributions [|30| almost 



coincides with the result of our model. As seen from the figure these ratios do not provide 
a test sensitive enough. 

An alternative idea to study the asymmetry was proposed more than a decade ago by 



Ito et al. [|T!|. They have suggested to analyze the logarithmic derivative of the rapidity 
distribution 



S(V^) = 4- In- d2fT 



, with r = Xix 2 , (27) 



dy dMdy 

where y = ln(x\/ x 2 ) /2 is the rapidity. This quantity also possesses the desired property 
of being independent of the .fT-factor. In terms of the quark distributions the slope can be 
expressed as 
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where 

4 Z_ N- _ Z N 

X(xi, x 2 ) = - {-^u{x 2 ) + -jd(x 2 )) + u(xt) [-ju{x 2 ) + -jd(x 2 ))} (29) 

+ ^ {d{xi)[-^d{x 2 ) + —u(x 2 )] + d(xi)[— d(x 2 ) + ^u(x 2 )]} 

+ n { s ( x i)^M + s(x 1 )s(x 2 )} . 
y 

We illustrate the effect of the u — d asymmetry on the slope of the rapidity distribution 
in Fig. 6. Here the solid lines are calculated using the asymmetric quark distribution of 
the recent MSR(A) quark parameterization |]30| and those obtained from our meson model 



|, |31|j3q| . The dashed lines are obtained by using symmetrized quark distributions given 
by Eq. (p3[) . The arrows show the effect of the symmetrization which decreases the slope. In 
this and all following calculations we have included corrections due to the Q 2 scale changing 
(Q 2 = sxix 2 ) by employing leading log Altarelli-Parisi QCD evolution. It has turned out 
that the resulting effects of the evolution are rather small. 

The rapidity slope (see Eq. |27|) is a quantity which is sensitive not only to the u — d 
asymmetry but also to valence quark distributions. In Fig. 7 we display the slope of the 
rapidity distribution calculated with different quark distributions. The solid line is the 
result of our meson model 1361. The dotted line is the result obtained with the Owens 



parameterization |40| of the quark distributions, the dashed line was obtained with the 
recent MRS (A) parameterization |3(| with u — d asymmetry and the dash-dotted line was 
obtained with MRS(S'q) 0] (symmetric) distribution. Fig. 7 clearly demonstrates that the 
asymmetry is not the only ingredient and a reasonable description of the experimental data 
|T2|j can be obtained with both flavour symmetric and asymmetric distributions. 



In obtaining both Eq.(p5|) and Eq. (|28|) we have neglected completely all nuclear effects 
like Fermi-motion, nuclear binding, excess pions or shadowing. Although they are predom- 
inantly flavour symmetric, it is obvious that they can modify the conclusions drawn based 
on the nuclear data. An information on nuclear effects can result only from the comparison 
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of the slopes for a nuclear target and for a deuteron one. It seems essential in the future 
to analyze more elementary processes, i.e. dilepton production in the proton-proton and 
proton-deuteron collisions. In the following we shall concentrate on those reactions as most 
reliable source of the information on the flavour asymmetry of the sea quarks. 
A quantity which can be extracted almost directly from the experiment is 

A DY (x u x 2 ) = , (30) 

GppKXli x 2) T &pn\ x l] x 2) 

which we will call Drell-Yan asymmetry. In formula (j3(|) o pp and a pn are the cross sections for 



the dilepton production in the proton-proton and proton-neutron scattering. The Drell-Yan 
asymmetry (|30"D can be expressed in terms of q and A introduced in Eq.(^j) 

_ [u(x 2 ) - d(x 2 )][3g(x 1 ) - 5/2A(xi)] - [4m(si) - d(x 1 )]A(x 2 ) , . 

dy{x 1 ,x 2 ) + rf(x2)][5 - (xi) _ 3 /2A(xi)] + [4w(xi) + d{x x )\2q{x 2 ) ' [ ' 

In the case of flavour symmetric sea (A = 0) it is natural to expect that A^y > since 
u > d. The sign of Apy can be, however, reversed by increasing the flavour asymmetry of 
the proton sea (A > 0). 

Two-dimensional maps of the Drell-Yan asymmetry as a function of Xi and x 2 are shown 
in the form of the contour plots in Fig. 8. The different maps were obtained with the Owens 
parameterization |^0| (left-upper corner), symmetric MRS(5 , g) @ (right -upper corner), the 
new MRS (A) 0] with the u — d asymmetry built in (left-lower corner) and the prediction of 



the meson cloud model [ 27|,3l|,p6 | (right-lower corner). The result obtained with the Owens 
(symmetric) parameterization and symmetric MRS(S'g) parameterization are quite similar. 
This clearly demonstrates that the asymmetry Apy is the desired quantity - insensitive to 
the valence quark distributions. It is also worth noting here that Apy is positive in the 
whole range of (xi,x 2 ). How the u — d asymmetry influences A^y is shown in two lower 
panels. It is very promising that A^y obtained with the asymmetric quark distributions 
(lower panels) differs considerably (please note the change of sign in the lower panels) from 
the result obtained with symmetric distribution (upper panels) and this should make an 
unambiguous verification of the flavour asymmetry of the sea quarks possible. It is not 
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random in our opinion that the result obtained within the meson cloud model is very similar 
to that obtained from the parameterization fitted to different experimental data. We stress 
in this context that A^y calculated in the meson cloud model is fairly insensitive to the quark 
distributions in the bare nucleons (baryons). It is primarily sensitive to the u — d asymmetry 
which is fully determined by the quark distributions in pions (mesons), taken here from the 
pion-nucleus Drell-Yan processes. We have assumed that the quark distributions in other 
mesons are related to those for the pion via SU(3) / symmetry. 

Following the suggestion of Ellis and Stirling, the NA51 Collaboration at CERN has 
measured recently the Ajyy asymmetry along the x\ = X2 diagonal ||. Due to low statistics 
only Aoy at low x = x\ = X2 was obtained. In Fig. 9 we show their experimental result (one 
experimental point) together with the results obtained with different quark distributions. 
The meaning of the lines here is the same as in Fig. 6. The result denoted as MCM, obtained 



within the meson cloud model [3j]|27 , |36| essentially without free parameters, nicely agrees 
with the experimental data point. In order to better understand the result and the relation 
to the u — d asymmetry let us express the cross sections in Eq. ( [30]) in terms of the quark 
distributions. Assuming proton-neutron isospin symmetry and taking X\ = X2 = x as for 
the NA51 experiment one gets in terms of the quark distributions in the proton 

5(« — d)(u — d) + 3(uu — dd) 
DY ~ 5{u + d){u + d) + 3{uu - dd) + k{ss + Acc)' 

Let us consider first the case u — d. For a crude estimation one may neglect sea-sea terms 
(important at small x only) and assume u va i(x) = 2d va i(x), which leads to A DY = 1/1 1 > 0. 
The same crude estimate in the case of asymmetric sea in conjuction with decomposition 
Eq.(|) yields 

= -19A + 6g 

-9A + 66g V ; 

This demonstates a strong sensitivity both ond-M asymmetry and on the absolute nor- 
malization of the sea. The lack of dependence on the valence quark distributions in the 
approximate Eq.(|55p suggests a weak dependence in the exact formula Eq. (|52"D . The nega- 
tive value obtained by NA51 experiment A^y = —0.09 ± 0.02 ± 0.025 automatically implies 
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d > u at least for the measured x = 0.18 (provided that the proton-neutron isospin symme- 
try violation is small(!)). The data point of the NA51 group is up to now the most direct 
evidence for the flavour asymmetry of the sea quarks, which is explicitly shown in Fig. 10 
where A DY has been translated into the ratio of u(x) /d(x) . The x dependence of the asym- 
metry is awaiting further experiments. It is expected that the new experiment planned at 



Fermilab [13] will be very useful in this respect and will provide the x dependence of the 
u — d asymmetry up to x = 0.4 and will shed new light on the microscopic structure of 
the nucleon. The meson-cloud model gives definite predictions for the asymmetry awaiting 
future experimental verification. 
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V. CONCLUSIONS 



The violation of the Gottfried Sum Rule observed by NMC 0|| together with nega- 
tive Drell-Yan asymmetry measured recently || by the NA51 group at CERN give a support 
to the conclusion that the SU (2) symmetry of the nucleon sea is violated. As discussed re- 
cently by Forte [IB] there are two possible kinds of symmetry violations, called SU(2)q and 
SU(2)i. The first one is simply connected with the asymmetry of light sea antiquarks u — d 
in the proton. The second is related to the violation of the proton-neutron isospin symmetry. 
There are no a priori reasons, except of customs of practitioners in deep-inelastic scattering, 
for either symmetry to be more fundamental. Both GSR violation and negative Drell-Yan 
asymmetry can in principle be explained by either SU(2)q (d > u) or/and SU(2)i (more 
abundant neutron than proton sea) symmetry violation. Some theoretical arguments fl6| 



suggest, however, that the violation of the SU(2)q symmetry seems to be more probable. 
While at present models explaining the u — d asymmetry have been constructed, no reliable 
models explaining the proton- neutron isospin symmetry violation exist. The proton- neutron 
symmetry violation can be expected on the basis of the bag models as due to the mass dif- 
ference of the u and d quarks as well as the corresponding di-quark states. At present a 
reasonable results can only be obtained for the valence quarks p5| , which, however, has no 
influence on the GSR violation and rather little effect for the Drell-Yan asymmetry at the 
experimentally measured x ~ 0.2. 

The old concept of the meson cloud in the nucleon gives a natural explanation of the u — d 
asymmetry. The essential parameters of the model - coupling constants - are well known 
from the low-energy physics. If the remaining parameters of the model (cut-off parameters 



of the vertex form factors) are fixed from the high-energy neutron and A ++ production [27 



then both the GSR violation and the Drell-Yan asymmetry can be well described. The same 
model gives also a good description of the neutron electric form factor |PJ] . In our model the 
virtual mesons influence the static properties of the nucleon (axial vector coupling constant 



electromagnetic radii etc.) 
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Our model has to be contrasted to the solution of Ball and Forte p5U46[| where mesons are 
produced radiatively via modified Altarelli-Parisi equations. Therefore their approach pre- 
dicts strong dependence of the Gottfried Sum Rule on the scale, in the range of intermediate 
Q 2 . In our approach the Gottfried Sum Rule is constant, at least in the leading logarithm 
approximation. It would be very important to test these two scenarios experimentally. A 
preliminary evaluation of the NMC data seems to support rather our model j|7| . 

We have discussed possibilities to identify the u—d asymmetry through the observation of 
the dilepton pairs in the hadronic collisions. The analysis of the present [i + pair creation 
data in proton-nucleus scattering is not conclusive. The rapidity slope, very sensitive to the 
flavour asymmetry, depends also on the valence quark distributions. The ratio of the cross 
section in proton-nucleus to that in proton-deuteron collision is compatible with symmetric 
sea quark distributions. However, the case of asymmetry concentrated at rather small x 
is not excluded. The meson cloud model gives results compatible with the E772 Fermilab 
experimental data [fTTfl . Elementary nucleon-nucleon Drell-Yan processes seem to be much 
better suited to study the d — u asymmetry. 

The presence of virtual mesons in the nucleon, especially pions, can explain the new 
result of the NA51 group at CERN for the Drell-Yan asymmetry. The Drell-Yan asymmetry 
is a quantity fairly insensitive to the valence quark distribution and very sensitive to the 
flavour asymmetry of the sea. The Drell-Yan asymmetries obtained with different valence 
quark distributions and symmetric sea are similar and positive. The meson cloud model 
predicts negative Apy, which is consistent with the only existing experimental data point 
||. Furthermore it gives definite prediction for the x dependence of the flavour asymmetry, 



awaiting experimental verification. The new experiment planned at Fermilab |T4[ will open 
such a possibility. 



One of us (A.S.) is very grateful to Peter Sutton for a very useful discussion on higher- 
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order QCD corrections. M.E. thanks the Humboldt Foundation for her stay in Germany. 
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FIGURES 

FIG. 1. The x-dependence of (x) — (x) . The solid line was obtained including the effect 
of virtual mesons. The data points at Q 2 = 4 (GeV/c) 2 are taken from Ref. ||. 

FIG. 2. The total sea quark distribution x(u(x) + d(x) + s(x)) obtained from our model. The 
solid line is the result of the full calculation, i.e. including the sea quarks in the bare nucleons 



(see Eq. 19). For comparison by the dashed line we show result obtained when neglecting the sea 
quark distribution in the bare baryons (mesonic contribution). The experimental data was taken 
from ||§. 

FIG. 3. The .ff-factors for the dilepton production in pp (solid line) and pn (dashed line) 
collisions for two different leading order quark distributions taken from [^0] and JXT| . 

FIG. 4. The cross section M 3 d 2 a/dxpdM for the production of the dilepton pairs in the 
proton-deuteron collisions. Shown is the fit of the K-factors for various quark distributions to the 



experimental data [29] 



FIG. 5. The Drell-Yan ratio defined by Eq.(|l3|) for iron and tungsten. The solid line is the 
result of our model, the dashed line the MSR(A) parametrization and the dash-dotted line the 
MSR(S'q) parameterization. 

FIG. 6. The effect of the u — d asymmetry on the slope of the rapidity distribution. The two 
solid lines were calculated with the asymmetric quark distribution of the recent MRS(A) quark 
parameterization [30] and that obtained from our model (MCM) |31. 27,36]. The dashed lines were 



obtained using the symmetrization procedure (see Eq. |23| ). The arrows show the effect of the 
symmetrization in both cases. 
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FIG. 7. The slope of the rapidity distribution. The solid line is the result of our model [36] 



The dotted line is the result obtained with the Owens parameterization |4C], the dashed line was 
obtained with the MRS(A) parameterization [|30| with d — u asymmetry and the dash-dotted line 
was obtained using the MRS(S ) (symmetric) distribution. 

FIG. 8. A two-dimensional map of the Drell-Yan asymmetry as a function of x\ and x^- 
Shown are results obtained with the Owens parameterization |2(| (left-upper corner), symmetric 
MSR(S ) (right -upper corner), new MRS (A) with the u — d asymmetry built in (left-lower corner) 
and prediction of the meson cloud model pl| , p^j36] | (right-lower corner). Note the change of the 
sign in the lower panels. 

FIG. 9. The Drell-Yan asymmetry along the x\ = X2 diagonal. The meaning of the lines here 
is the same as in Fig. 6. The experimental data point is taken from the recent result of the NA51 
Collaboration at CERN 0. 



FIG. 10. The u(x)/d(x) ratio as obtained from the meson cloud model (solid line) |3(| and 
the MSR(A) parameterization (dashed line) compared with the experimental result of the NA51 
collaboration 101. 
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